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We study a parametrically modulated oscillator that has two stable states of vibrations at half the 
modulation frequency luf- Fluctuations of the oscillator lead to interstate switching. A compara- 
tively weak additional field can strongly affect the switching rates, because it changes the switching 
activation energies. The change is linear in the field amplitude. When the additional field frequency 
LOd is luf/2, the field makes the populations of the vibrational states different thus lifting the states 
symmetry. If ujd differs from uf/2, the field modulates the state populations at the difference fre- 
quency, leading to fluctuation-mediated wave mixing. For an underdamped oscillator, the change 
of the activation energy displays characteristic resonant peaks as a function of frequency. 



PACS numbers: 05.40.-a, 05.70.Ln, 74.50.+r, 02.50.-r 



I. INTRODUCTION 

A parametrically modulated oscillator is one of the 
simplest physical systems that display spontaneous 
breaking of time-translation symmetry. When the mod- 
ulation is sufficiently strong, the oscillator has two states 
of vibrations at half the modulation frequency, the pe- 
riod two states pj. They are identical except for the 
phase shift by n, but for each of them the symmetry with 
respect to time translation by the modulation period is 
broken. Fluctuations of the oscillator lead to switchings 
between the states. Switching rates are equal by sym- 
metry, for stationary fluctuations. The switchings ulti- 
mately make the state populations equal, thus restoring 
the full time-translation symmetry. Experimental studies 
of fluctuation-induced switchings in classical parametri- 
cally modulated systems were done for trapped electrons 
Q , optically trapped atoms 0, 3 , and microelectrome- 
chanical systems |5j . The obtained switching rates are in 
good agreement with the theory 

The degeneracy of period two states can be lifted if, 
in addition to parametric modulation at frequency wp, a 
system is driven at frequency lof/2. In the frame oscillat- 
ing at frequency ujf/2 the period two states and the addi- 
tional field look static. The system reminds an Ising fer- 
romagnet, with the period two states and the additional 
field playing the roles of spin orientations and an exter- 
nal magnetic field, respectively. One can expect that the 
role of the direction of the magnetic field is played by the 
phase of the additional field counted off from the phase 
of one of the period two states. Depending on the field 
phase, one or the other state should be predominantly 
occupied. 

In this paper we study resonant symmetry lifting in a 
parametric oscillator, which occurs where the frequency 
of the additional field is close to the oscillator eigenfre- 
quency. The field makes the rates of switching between 
the period two states, Wn and W21, different from each 
other. In turn, this leads to the difference of the sta- 
tionary state populations. This difference may become 
large even for a comparatively weak field, as can be sur- 
mised from the analogy with the problem of a fcrromag- 



net. There the change of the state populations becomes 
large once the energy difference of the states due to an 
external magnetic field exceeds kT, which happens al- 
ready for weak fields where this difference itself is small 
compared to the internal energy. 

In contrast to the case of a ferromagnet, the energy 
of a parametrically modulated oscillator is not conserved 
and its stationary distribution is not of the Boltzmann 
form. However if the fluctuation intensity is small, the 
dependence of the switching rates on this intensity is of- 
ten of the activation type, see Refs.0,0and papers cited 
therein. This applies not only to a classical, but also to 
a quantum oscillator, where switching is due to quantum 
fluctuations Q . An additional field changes the effective 
switching activation energies, and once this change ex- 
ceeds the fluctuation intensity, the overall change of the 
switching rates becomes large. The effect is particularly 
strong if the field is resonant. 

In what follows we develop a theory of the switching 
rates W nm for a classical oscillator. We find the depen- 
dence of the switching activation energies on the ampli- 
tude, phase, and frequency of the additional field and on 
the oscillator parameters. We first study the symmetry 
lifting by the field at frequency u>d = u>p/2. Of partic- 
ular interest here is the vicinity of the bifurcation point 
where the period two states merge, as in this range the 
rates W nm are comparatively large and easy to control. 

We are also interested in the situation where the ad- 
ditional field frequency ujd is not exactly equal to ujp/2. 
Here, the field-induced modulation of the switching prob- 
abilities causes oscillations of the state populations at fre- 
quency \ujd — luf/2\. Such oscillations, superimposed on 
the oscillator vibrations at frequency luf/2 in period two 
states, lead to vibrations at frequency \u>f — u>d\, i.e., to 
a strong effective fluctuation-induced three- wave mixing. 

The amplitude of the population oscillations becomes 
small when \u>d — LOpj2\ largely exceeds the switching 
rates. Still the rates themselves may be significantly 
changed by the additional field. Of primary interest in 
this case are the rates W nm averaged over the period 
47r/|2wd — u>f\ and their dependence on the additional 
field amplitude Ad and frequency Ud- One may expect 
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that, as in the case of equilibrium systems [ToL ITTl Il2^. 

the rate change is quadratic in A d for small Ad and 
corresponds to effective heating of the system by the 
field. For underdamped equilibrium systems this heat- 
ing can be resonant, as seen for modulated Josephson 
junctions However, for somewhat stronger fields 

the change of the logarithm of the switching rate should 
become linear in Ad |l5l ITr| . This happens when the 
properly scaled field amplitude exceeds temperature. 

We show that, for a parametric oscillator, the acti- 
vation energies are indeed linear in the additional field 
amplitude Ad once it is not too small (but is also not 
too large) . When the oscillator is underdamped, the fac- 
tor multiplying Ad displays a characteristic resonant fre- 
quency dependence. We develop a technique that allows 
us to find this dependence in an explicit form, for the 
considered nonequilibrium system. The asymptotic ana- 
lytical results are compared with the results of numerical 
calculations of the activation energies. 

In Sec. II we discuss the Langevin equation for a para- 
metrically modulated nonlinear oscillator in the rotating 
frame and give the general expression for the probability 
of switching between coexisting stable vibrational states. 
In Sec. Ill we obtain a general expression for the cor- 
rection to the activation energy of switching and show 
that it is linear in Ad- We study resonant symmetry 
lifting of the switching rates. In Sec. IV we investigate 
low-frequency oscillations of the state populations and 
fluctuations-mediated resonant wave mixing. In Sec. V 
an expression for the period-averaged switching rate is 
given. In Sec. VI we consider symmetry lifting and the 
frequency dependence of the activation energy close to 
the bifurcation point where the period two states merge. 
In Sec. VII we study the case of weak damping. We 
show that the change of the activation energy may dis- 
play characteristic asymmetric resonant peaks as a func- 
tion of the additional field frequency. Sec. VIII contains 
concluding remarks. 



II. LANGEVIN EQUATION AND THE 
SWITCHING RATES 

We will study switching between period two states of 

a nonlinear oscillator, which is parametrically modulated 

by a force Fcos(uj F t) and is additionally driven by a 

comparatively weak field AdCOs(u> d t + 4>d) at frequency 

bJd ~ uj f /2. The Hamiltonian of the oscillator is a sum 

of the term that describes the motion without the extra 

field and the term proportional to the field, H osc = Hosl + 
H (d) 

H { £> c = -pl + \ql{ul + F C os{u F t))+ l - lq t, 
H^l = -q A d cos{u d t + <f> d ) (1) 

(<?o and po are the coordinate and momentum of the os- 
cillator). We will assume that the modulation frequency 



ujp is close to twice the frequency of small amplitude vi- 
brations Wo, and that the driving force F is not too large 
so that the oscillator nonlinearity remains small, 

|wp - 2w | , \uj d - uj \ <C uj , (2) 
F<£u, 2 , | 7 |<g 2 )«t4 

In what follows for concreteness we set 7 > 0. 

Following the standard procedure pj , we change to the 
rotating frame and introduce the dimensionless canonical 
coordinate Q and momentum P, 

q (t) = C[Pcos{uj F t/2) - Qsu\(uj F t/2)] , (3) 
Po (t) = -C^-[Psm(u F t/2) + Qcos(uj F t/2)}, 

where C = (2F/3"/) 1 / 2 . In these variables the Hamil- 
tonian becomes equal to H osc = (F 2 /6j)g(Q, P), with 
g = gl°)+gW(T), 

g {d) {T) = -a d [Pcos(v d T + <f> d ) 

+Q sin (v d T + <t> d )} . (4) 

Here we introduced dimensionless time r and dimension- 
less parameters [i and v d . These parameters characterize, 
respectively, the detuning of the modulation frequency 
from the oscillator cigenfrcquency and the detuning of the 
weak-field frequency from uj f /2, i.e., an effective "beat 
frequency" with the subharmonic of the strong field, 

fi = lo f (uj f - 2u )/F, v d = u F (2u> d - uj f )/F, 

t = tF/2uj F . (5) 

The parameter a d = Ad{&^ / F 3 ) 1 / 2 is the dimensionless 
amplitude of the additional driving field. In obtaining 
Eq. we used the rotating wave approximation and 
disregarded fast oscillating terms cx exp(±inLu F t) , n > 1. 
In the quantum formulation, the eigenvalues of </°) (Q, P) 
give the scaled quasienergy of the system , and in what 
follows for brevity we call g quasienergy. 

We will assume that the interaction with a bath that 
leads to dissipation of the oscillator is sufficiently weak, 
so that the oscillator is underdamped. Then under fairly 
general assumptions [l7j j. in the rotating frame dissipa- 
tion is described by an instantaneous friction force (no 
retardation). Also the noise spectrum is generally prac- 
tically flat in a comparatively narrow frequency range of 
width > F/lu f centered at o>o; this is the most interesting 
range, since the oscillator filters out noise at frequencies 
far from this range. Therefore with respect to the slow 
time t the noise can be assumed white. The oscillator 
motion is described by the Langevin equation, which can 
be conveniently written in a vector form as 

q=^=K + f(r), K = KW+K( d ), (6) 
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with 

K<°> = -C _1 q + eV«7 (0) , K (d >(r) = eV. 9 (d) (r). (7) 

Here all vectors have two components, q = (Q,P), 
K = (Kq, Kp), and V = (8q, dp), while e is the permu- 
tation tensor, eqq = epp = 0, egp = — epQ = 1. The 
parameter in Eq. (JSJ gives the oscillator friction co- 
efficient in the units of F/2top. We use for ( and fi the 
same notations as in Refs.|a,0- Because of the extra field 
oc ad, the function K^ d ' explicitly depends on time. 

The function f (r) is a random force. Its two compo- 
nents are inde pen dent white Gaussian noises with the 
same intensity |l7| . 

(/q(t)/q(0)) = (fp(r)fp(O)) = 2DS(r - t'). 

The noise intensity D is the smallest parameter of the 
theory. If both the friction force and the noise come 
from coupling to a thermal reservoir at temperature T, 
we have D = 67/cT JF^oj f (the parameter D corresponds 
to L>C~ 2 /2 in Ref.@). . 

A. Oscillator dynamics in the absence of noise 

In the absence of noise and the symmetry-breaking 
field oc Ad, in the range 

^ 2 ) = T( i-r 2 ) 1/2 (8) 

the parametrically modulated oscillator has two stable 
period- two states qf \ and an unstable state q^ . These 
states are the stationary solutions of equation K' ' = 0. 
They merge for fj, = f/g . The stable states 1, 2 are 

inversely symmetrical, q^ = — . For concreteness 
we choose 

= -QP > 0. 

The vibration amplitude in the unstable state is zero, 
q^ = 0. For \x > fig the state q = becomes stable 
and there additionally emerge two unstable period two 
states 0. 

We will be interested in a comparatively weak 
symmetry-breaking field. Respectively, we will assume 
that the reduced field amplitude ad is small, so that the 
field does not lead to new stable states. It just makes the 
stationary states periodic, for vj, ^ 0, or shifts them, for 
Vd = 0. The correspondingly modified states are given by 
the periodic solutions of equation q = K or by equation 
K = 0. In the laboratory frame the periodic stable states 
qi.2(i~) correspond to oscillator vibrations at frequency 
uip/2 weakly modulated at frequency \u>p — 2u)d\/2. They 
have spectral components at ujp/2,ujd, and the "mirror" 
frequency \u>p — u>d\- 

We will limit ourselves to the analysis of switching in 
the parameter range J5J). In this range escape from a 



period two state leads to switching to a different period 

(2) 

two state. For u > yr B escape may result in a transition 
to the zero-amplitude state (from which the system may 
also escape to one of the period two states) . The results 
of the paper immediately extend to this range, but this 
extension will not be discussed. Therefore we use the 
terms "escape" and "switching" intermittently 

B. Switching rates: General formulation 

The noise f(i) leads to fluctuations about the stable 
states and to interstate transitions. When the noise is 
weak, fluctuations have small amplitude on average. In- 
terstate transitions require large outbursts of noise and 
therefore occur infrequently For Gaussian (5-correlatcd 
noise the probability of a transition from nth to mth sta- 
ble period two state has activation dependence on the 
noise intensity D and is given by the expression 0, @ 

W nm = C w eM-Rn/D), (9) 

f°° 1 
R n = min / dri(q, q; t) , L=-(q-K) 2 . 

J —00 

The quantity R n is the activation energy of a transition. 
It is given by the solution of a variational problem. The 
minimum in Eq. 10 for R n is taken with respect to tra- 
jectories q(r) that start for r — > —00 at the initially oc- 
cupied stable state q ra (T) and asymptotically approach 
qt,(r) for r — > 00. The optimal trajectory q n0 pt(*) that 
minimizes R n is called the most probable escape path 
(MPEP). The system is most likely to follow this trajec- 
tory in a large fluctuation that leads to escape. 

The prefactor in the transition probability is Cw oc 
F/ujp in unsealed time t. It weakly depends on the noise 
intensity and will not be discussed in what follows. 

The variational problem © for R n can be associated 
with the problem of dynamics of an auxiliary system 
with Lagrangian L and coordinates q = (Q, P), and with 
Hamiltonian 

h = h<-Q+hW, ff (0) = p 2 + pKW, (10) 

In terms of this auxiliary Hamiltonian system, the MPEP 
of the original dissipative system corresponds to a het- 
eroclinic trajectory that goes from the periodic (or sta- 
tionary, for Vd = 0) state (q n (T),P = 0) to the periodic 
state (qh(r), p = 0). 

III. PERTURBATION THEORY FOR THE 
ACTIVATION ENERGY 

When the additional field oc ad is weak, the activation 
energy R n is close to its value R^ for aa = 0, which is 
the same for the both states n = 1,2. Even though the 
correction to R^ is small compared to R^ , it may sig- 
nificantly exceed the noise intensity D. Then the overall 
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change of the transition probability W nm will be expo- 
nentially large. 

The correction to the activation energy was studied 
earlier for thermal equilibrium systems additionally mod- 
ulated by a comparatively weak periodic field 0, 
As mentioned in the Introduction, it was found that the 
correction to R n is linear in the field amplitude. The 
factor multiplying the amplitude gives the logarithm of 
the transition probability and therefore was called the 
logarithmic susceptibility (LS). We show now that the 
correction to the activation energy for a parametrically 
modulated oscillator is also linear in the amplitude of 
the additional field, SR n = R n — R(°> oc a<j> an d nn d the 
proportionality coefficient, that is the LS. 

Because of the modulation at frequency ojf, the oscil- 
lator is far away from thermal equilibrium even without 
the field oc ad- This leads to a significant difference of 
the problem of switching from that for equilibrium sys- 
tems. Its physical origin is the lack of time reversibil- 
ity in noncquilibrium systems. As a result the auxiliary 
Hamiltonian system described by the Hamiltonian 
()10JI . which gives optimal fluctuational trajectories and 
the MPEP of the original dissipative system, is nonintc- 
grable [is[ . 

The Hamiltonian dynamics described by the full 
Hamiltonian H l)l()|l. which includes the additional time- 
dependent field cx ad, is also nonintegrable. However, for 
small ad the heteroclinic Hamiltonian trajectory q rl0 pt(''") 
that gives the MPEP may remain close to the unper- 
turbed heteroclinic trajectory q„o pt (r) for ad = 0. Then, 
as in the case of modulated equilibrium systems |l5l Il6j| , 
the lowest-order correction in ad to the variational func- 
tional R n © can be found by calculating the perturbing 
term in the Lagrangian along q^l pt (r). 

The trajec tory qi°o P t( T ) is a real-time counterpart of 
an instanton |l9| and is often called a real-time instanton. 
It goes from r — > — oo to r — > oo. As in the case of 
standard instantons, q^ Q pt (r) looks like a pulse. It is 
large only for a time of the order of the relaxation time 
of the system in the absence of noise. The position of the 
center of the pulse on the time axis r c is arbitrary. 

Periodic modulation lifts the degeneracy with respect 
to t c . It synchronizes switching events. It is this synchro- 
nization that leads to the correction to R n being linear 
in ad for a nonzero frequency detuning \vd\- The syn- 
chronization corresponds to calculating the field-induced 
correction 8R n using as a zeroth order approximation the 
unperturbed MPEP for the nth state qi°opt( r — T c) and 
adjusting r c in such a way that the overall functional 
Rn = i?*- ' + 5R n be minimal. That is, r c is found from 
the condition of the minimum of the function SR n (T c ). 



SRn 
SRn(T c ) 

Xu{t) 



min 5R„(t c ), 

/OO 
-OO 



(11) 



Clearly, r c is determined modulo the dimcnsionless mod- 
ulation period 2n/\i'd\- Generically, there is one MPEP 
per period that provides the absolute minimum to R n . 

Equation (|11|) shows that the correction to the acti- 
vation energy of escape is indeed linear in ad- The co- 
efficient at ad is determined by the Fourier components 
Xn(±v d ) of the function x(t), 



r Xn(r). 



(12) 



The function Xn ( r ) an d its Fourier transform Xn (y) give 
the LS in the time- and frequency representation. They 
are of central interest for the studies of symmetry lift- 
ing and switching rate modulation. Interestingly, the 
functions Xn ( r )i Xn(y) have two components each, even 
though there is only one additional force that drives the 
oscillator. This leads to important consequences dis- 
cussed below and is in contrast to the case of a modulated 
equilibrium system. 

We note that the condition d[SR n (T c )]/ 'dr c = 0, 
Eq. (jl 1|> , corresponds to the condition that the MePnikov 
function M.(t c ) defined for our Hamiltonian system with 
two degrees of freedom as 



M{t c ) 



[h<® (T-r c ),H»(r-r c ;r)},(13) 



dr{H 



qi°) pt (r)-K(°) (qi°) pt (r 



be equal to zero. Here, {A, B} is the Poisson bracket 
with respect to the dynamical variables of the auxiliary 
system q, p. The Poisson bracket is evaluated along the 
unperturbed trajectory qi°o pt (r - t c ), pi°o pt (r - r c ). Re- 
spectively, the functions q, p in H^°\H^ are evaluated 
at time r — r c , which is indicated by the first argument 
in 

#(0) 5 #(<Q. t he second argument in flw indicates the 
explicit time dependence of H^ d > , Eq. I|l()(l . due to the 
field oc ad- For systems with one degree of freedom, the 
condition M(t c ) = shows that r c for the unperturbed 
trajectory is chosen in such a way that this trajectory is 
close to the trajectory q„ op t(i" - r c ), p nop t(T - r c ) [2(|. 
For our system the condition M(t c ) = is necessary for 
the applicability of perturbation theory. The correspond- 
ing analysis will be provided elsewhere. 



A. Symmetry lifting by a field at subharmonic 
frequency 

The analysis of the switching rates should be done 
somewhat differently in the case where the frequencies 
of the additional field and the parametrically modulating 
field satisfy the condition ujd = luf/2- In this case Ud = 
and the force K*-^ is independent of time. Therefore the 
function SR n (T c ) is independent of r c and minimizing it 
over t c is irrelevant. However, a perturbation theory in 
ad can still be developed. The first order correction to 
the activation energy SR n = SR n cs (4>d) can be obtained 
by evaluating the term oc ad in the Lagrangian L @ 
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along the unperturbed path q^ Q pt (r). It has a simple 
explicit form 

6R™ s ((t) d ) = a d [x, l q(O)cos0 (i - x„ P (O)sin0 d ] , (14) 

where the subscripts Q,P enumerate the components of 
X [we have also used Eqs. (J3J, J7J) for K^]. Both compo- 
nents of Xn(0) contribute to the change of the activation 
energy for subharmonic driving. 

An important feature to emphasize is that the function 
5R I 1 ^ s (4>d) is not the limit of min 8R n (T c ) for v& — > 0. 
Physically this is because Eq. (|11J) gives the change of 
the logarithm of the escape rate W nm averaged over the 
period ir/\2uJd — u>p\ (or 2ir/\i'd\, in dimensionlcss time 
r). Such averaging is meaningful as long as the frequency 
detuning \2uj d — ujp\ >• W nm . In the opposite limit the 
occupation of the state n changes significantly over the 
period ir/\2wd — ujp\; this change is not characterized by 
the period-averaged W nrn . 

An important property of the static LS Xn (0) is that 
it has opposite signs for the states 1 and 2. Indeed, these 
states are inversely symmetric, qj = — q^ ' , and since 
the vector K^ ) is also antisymmetric, K^ ) — > — K^ ) 
for q — > q, it is clear that c^l t {T) = — ^optC 7 "); an d 
therefore Xi(0) = -%2(0). Then from Eqs. ©, JE} we 
have for the switching rates W\ 2 and W 2 \ 

W 12 /W {0) = W w /W 21 = exp [-8Rl es (<f> d )/D] , (15) 

where = = W^i i s the switching rate for 

a d = 0, WW oc exp (-R.W/D). Eq. (O applies for 
arbitrary ad/D. The corrections to the prcfactor cx ad 
(but not ad/D) have been disregarded. 

In the parameter range fig < fi < fig where the 
only stable states of the system are the period two states 
[cf. Eq. ©], the ratio of the stationary populations of the 
states wi and w 2 is the inverse of the escape probabilities 
ratio. Therefore from Eq. I|15|) 

Wl /W2 = W 21 /W 12 = exp^rOdV-D] ■ 



It is seen from this expression that even a compara- 
tively weak symmetry-lifting field, where l^i? 1 ^! -C R^°\ 
can lead to a significant change of the state populations. 
This happens when |<5i?i cs ^ D. The ratio of the state 
populations is determined by the phase of the field (fid- 
As mentioned in the Introduction, there is a similarity 
with magnetic field induced symmetry breaking for an 
Ising spin, with <f>d playing the role of the orientation 
of the magnetic field. We note that both the absolute 
value and the sign of the coefficients Xn(0) depend also 
on the oscillator parameters for ad = 0. Therefore one 
can change the population ratio not only by varying the 
phase and amplitude of the additional field, but also by 
varying these parameters, for example the amplitude F 
of the paramctrically modulating field. 



IV. LOW-FREQUENCY MODULATION OF 
STATE POPULATIONS 



The LS Xn(y) displays frequency dispersion for v of 
the order of the inverse dimensionless relaxation time of 
the oscillator or higher. If the additional field is so closely 
tuned to the subharmonic frequency that \2ujd — ojp\ < 
W(°\ the period averaging of W nm implied in Eq. i|ll|) 
becomes inapplicable, as explained above. One should 
rather think of the instantaneous values of the switch- 
ing rates W nm (T), which are given by Eq. (|15f) with 
4><i — > VdJ + 4>d(0). Slow time-dependent modulation of 
the switching rates leads to modulation of the state pop- 
ulations, 



— — = -Wi2{T)Wi + W2l(T)W2; 

ar 



W\ + w 2 = 1. (16) 



Such modulation has attracted much attention in the 
context of stochastic resonance j^, . In contrast to a 
particle in a slowly modulated double-well static poten- 
tial that has been most frequently studied in stochastic 
resonance, here the stable states are fast oscillating, and 
strong modulation of their populations is induced by a 
high-frequency driving field. In this sense there is sim- 
ilarity with stochastic resonance in a resonantly driven 
oscillator with coexisting period one states (see Ref. l2lj) 
that was recently observed in experiment |23j . In con- 
trast to a resonantly driven oscillator, for a paramctri- 
cally modulated oscillator the populations of period two 
states are equal in the absence of extra field for any pa- 
rameter values. 

Oscillations of the state populations at frequency \vj\ 
(\u>d — u>p/2\ in the laboratory frame) have a compara- 
tively large amplitude oc ad/D, for small noise intensity. 
They lead to vibrations of the oscillator in the laboratory 
frame at frequencies LOd and \ujf — LOd\, with the average 
coordinate being 

(q (t)) « |2F/3 7 | 1/2 KW-^(t)] 



x \p( 0) cos(w F </2) - Q { ° ] sinw F t/2 



The vibration amplitude is much larger than the am- 
plitude of noise-free vibrations about attractors. This 
indicates strong fluctuation-induced effective three-wave 
mixing. Moreover, because oscillations of w\, 2 {t) are 
nonsinusoidal for large ad/D, there also occurs multiple- 
wave mixing. Further analysis of this effect is beyond the 
scope of the present paper. 



V. PERIOD-AVERAGED SWITCHING RATE 

For the difference frequency \2ujd - lu f \ > Wi^ the 
major effect of the additional field is the change of 
the switching rate averaged over the dimensionless pe- 
riod 27r/|z/(j|. It is given by the change of the activa- 
tion energy 8R n . From Eq. Ijllfl it follows that 5R n 
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is independent of time. From the symmetry relations 
q^ pt (r) = -q^ pt (r) and KW(r) = -K.W (r + tt^ 1 )) 
it follows that the minimization over r c in Eq. leads 
to SRi = (5i?2 ■ It is straightforward to show that 



5Ri = SR 2 = -a<ixic(vd), 

xicM = [[xiMI 2 -i*iMexiM] 1/2 . 



(17) 



From Eq. |17jl. the change of the activation energy is 
fully determined by the LS at the scaled frequency dif- 
ference I'd- The function xic is nonnegative. It displays 
a characteristic dependence on the oscillator parameters. 
As we show, it may have resonant peaks in the regime of 
small damping, £ 3> 1. 



VI. SCALING BEHAVIOR NEAR THE 
BIFURCATION POINT 

A. Symmetry lifting 

The dynamics of the oscillator is simplified near the 
bifurcation point /j^ where the period two states merge 
together (a supercritical pitchfork bifurcation [13] )■ Here, 
motion is controlled by one slow variable (soft mode) 
Q' = Qcos/3 + Psin/3, where (5 = \{-k — arcsinC" 1 ). 
From Eq. ©, to the lowest order in the distance to the 
bifurcation point <r\ = jx — jA^ the Langevin equation for 
this variable has the form 

Q' = -d Q ,u + f{T), u = ^ v cq' 2 

-^b } C 3 Q' 4 + Q'a d cos(^r + 4>d + 0), (18) 

where /'(r) is white noise of intensity 2D. 

The potential U in the absence of additional field has a 
familiar form of a quartic parabola (note that /z^ < 0). 
The values of the slow variable at the period two states 

correspond to the minima of U, Q12 = i ? / 1 ^ 2 C~ 1 j 
whereas the unstable zero-amplitude state is at the lo- 
cal maximum of U at Q'^ = 0. The activation energy 
= Ma]?7 2 /4C is just the height of the potential bar- 
rier AU ]g; the probability distribution near the pitch- 
fork bifurcation point was discussed in Refs. [3 EH- 

For an additional field at exact subharmonic frequency, 
the change of the activation energy is simply the change 
of the barrier height. With the chosen convention that 
> wc obtain from Eq. (|18|) 



,(0) 



5R[ CS = -SR[ C 



V /2 C" 1 adCos(^ + ( 9). (19) 



This expression shows that the correction to the activa- 
tion energy decreases as the system approaches the bi- 
furcation point, i.e., r\ = fi — /Lt^ decreases. However, 
the decrease of SR\ CS is much slower than the decrease of 
R(°>, Therefore the relative correction to the activation 
energy sharply increases as fi — > . It is important 



that the sign of 5R\ es , that shows which of the states 1 
and 2 is predominantly occupied, depends on (f>d — (3, that 
is, not only on the relative phase of the additional field, 
but also on the scaled relaxation parameter £ that de- 
termines the value of (3. Therefore by varying £ one can 
control which of the states is predominantly occupied. 



B. High-frequency modulation 

We now consider the case where the additional field is 
detuned from the subharmonic frequency, \2ujd — uip\ 3> 
■ Since the motion near the bifurcation point is con- 
trolled by one variable Q', there is no phase shift between 
the components XiQ^d) and Xip{ v d)- The change of 
the activation energy is 5R\ = SR2 = —cidXic{vd) = 
~ a d\Xi{v d )\- The problem as a whole coincides with 
that for a periodically modulated overdamped equilib- 
rium particle. The LS for an overdamped particle in a 
quartic potential is already known ||. In the present 
case we have 

xicM = 7r-v»c-V /a 
r(d- 



((l-*V)/2)r(l+*V/2) , (20) 



where v' = u/\fj,g |rjC an d T( z ) is t ne Gamma function. 
The function X\ c {y) is proportional to a smaller power 
of the distance to the bifurcation point 77 than R^>. 
This shows that the correction to the period-averaged 
switching rate becomes relatively larger as the system 
approaches the bifurcation point. As seen from Eq. 12011 . 
for small 77 the LS Xiciy) has a peak at the frequency 
detuning v = and monotonically decreases with in- 
creasing \v\. The typical width of the peak of X\ c {v) is 
C|Mb l^- It decreases as /i approaches /U^ . 



VII. THE WEAK DAMPING LIMIT 

The LS can be analyzed also in the limit of weak damp- 
ing, ( > 1. Here we consider damping that is weak in 
the rotating frame. This means that not only is the oscil- 
lator decay slow compared to its frequency ujq, but also 
compared to the much smaller frequency F/uip- If there 
were no damping and noise, the motion of the oscilla- 
tor in the rotating frame would be vibrations with given 
quasienergy g, which are described by equation q = eVg. 
Damping causes the quasienergy to decrease towards its 
value in one of the stable states whereas noise leads to 
quasienergy diffusion away from these states. On the 
MPEP 

q« opt (t) the quasienergy increases from its value 
g n _in the stable state n to its value at the saddle point 
9b 



Q (the quasienergy G in Ref. differs from </°) (Q, P) 



in sign). As a result, q n opt(i) is a spiral. 
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A. Symmetry lifting 

We will start the analysis with the case of the addi- 
tional field at exact subharmonic frequency, ujd = ujp/2. 
In this case g(Q, P) is independent of time. The general 
expression for the activation energy of switching in the 
limit of weak dam ping , which is not limited to small ad, 
has the form lEHlIjl 



N n (g) = ^ 



9b , M n (g) 

9n d9 W' 



M n {g) 



dQdP, 



A n (g) 



V 2 g dQdP. 



A n (g) 



(21) 



Integration with respect to Q, P in functions 
M n (g), N n (g) for given g is done over the area A n (g) 
of phase plane (Q,P), which is limited by the phase 
trajectory g(Q, P) = g that lies within the basin of 
attraction to the stable state q„. 

For comparatively small a d the integrals in Eq. I|21|) 
can be calculated by perturbation theory. To first order 
in ad we obtain the following expression for the correction 



5R™ S = 5R™ s ((f>d) to the activation energy for ^ 



< 



where 
X 1 = 

SM 1 -- 

5N X -- 



SRT^d) 



dg 



\d? d ) = X 1 C 1 a d sin^ (22) 



Co) a AT(0) 
Qi 



6Mx tt^SN! 



AT(o) 



2 / QdQ/\d P g^\ 



(m + i) 1/2 



(23) 



QdQ V 2 g^/\d P g 



(0)| 



Here and are the values of M(g) and N{g) 

for ad = 0, whereas 5M and 5N are the field- induced 
corrections to M and N divided by a^sin^; g^( n = 
-(fj, + l) 2 /4. The limits Q lmin (g), Qi m ax(.g) of the inte- 
grals over Q are given by equation g^°\Q,0) = g, with 
< Qimin < Qimax- The arguments of the integrals 
over Q in the expressions for SMi , 5Ni are calculated for 
P given by equation g^(Q, P) = g with P > 0. 

It follows from Eq. 11221) that, for small damping, the 
dependence of the correction to the activation energy 
on the phase of the additional field has a simple form 
8R\ CS oc sm<?V Both 5R\ CS and are oc C"\ they de- 
crease with the decreasing scaled friction coefficient 
The function X\ depends on one parameter, the relative 
frequency detuning of the strong field [x. This depen- 
dence is shown in Fig. ^ 

It is seen from Fig. that the change of the activation 
energy is nonmonotonic as a function of \x. For /x close to 
the bifurcation value jjL^ = — 1 we have X\ w (fi+ l) 1 / 2 . 
This shows that the weak-damping expression for SR\ CS , 




FIG. 1: Symmetry lifting in the weak-damping limit, f~ <C 
1. Solid line: the scaling factor Xi, Eq. 123L in the correction 
to the activation energy SR{ CS = XiQ~ a^sinc^. The dots 
show the scaled LS CXf (0)- It gives the correction oc sin^d 
to 5R\ es and is obtained by calculating the MPEPs for ( = 5. 



Eqs. 1221) . I|23|) . smoothly goes over into expression 119|) 
obtained near the bifurcation point in the opposite limit 
of overdamped motion; note that in Eq. (|19|l j3 k, tt/2 for 

r 1 < 1. 

The analytical results are compared in Fig. ^ with the 
numerical results obtained by directly solving the Hamil- 
tonian equations of motion for the MPEP and calculat- 
ing the LS from Eqs. (fTT|) . 1(12)1, Since in the limit of 
small damping 5R\ CS oc sincA^, one expects from Eq. i|14|> 
that Xiq(0) should be small and X\ « Cxip(0)- Close 
to the bifurcation point fj,^ the numerical results agree 
with the asymptotic theory already for moderately small 
damping, = 0.2. 

We found numerically that the component Xiq(0) in- 

(2) 

creases with \i. Close to the second bifurcation point n B 
it becomes of the same order of magnitude as Xip(0), for 
the chosen £. We note that the MPEP is a fast oscillating 
function of time, for weak damping. When Xi(0) is cal- 
culated, the oscillations largely compensate each other. 
The numerical value of Xi(0) is therefore very sensitive 
to numerical errors in the MPEP. With increasing /z the 
frequency of oscillations of the MPEP increases and so 
does the sensitivity. Therefore we present numerical re- 
sults only for a limited range of fi. 



B. Resonant peaks of the logarithmic susceptibility 

Oscillations of the MPEP may be expected to lead 
to peaks of the LS as function of frequency for the di- 
mensionless frequency detuning \v\ » Since the 

period-averaged corrections to the activation energy for 
the states 1 and 2, SR\ and <5i?2, are equal, we will con- 
sider the MPEP for the state 1 and will drop the sub- 
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script 1. The analysis of the LS peaks can be done by 
extending to the MPEP of a parametrically modulated 
oscillator the approach developed in Ref . [T3| ■ 

As a first step it is convenient to change variables from 
Q,P to quasienergy-angle variables g,ip. This is accom- 
plished by seeking the optimal path of the unperturbed 
system, ad = 0, in the form 

qioptW = (^(ffo P t,V'opt),y(5opt, V'opt))- 

Here x , y are the coordinate and momentum of vibrations 
with given unperturbed quasienergy, g(°\x,y) = g; ip is 
the vibration phase: x and y are 27r-periodic in tJj. We 
denote the vibration frequency by u)(g). The equations 
for x , y are of the form 



v{g)d,px = d y g(°\ oj(g)d i ,y = -d x g 



(0) 



(24) 



Functions x{g,ip) and y(g,ip) can be expressed in terms 
of the Jacobi elliptic functions [|| . 

On the optimal path the quasienergy g = goptfV) is 
a function of time, with g opt , \tp op t - uj(g op t)\ oc Ca- 
using the explicit form of equations of motion for the 
Lagrangian one can show that, to first order in £ , 
on the optimal path 



(dx/dT) opt 
(dy/dT) opt 



K^(x,y) + C 1 F(g)d^y, 
K^(x,y)-C 1 F(g)d^x. (25) 

Here, the force K^, with components K x , Ky , is de- 
fined by Eq. Q in which Q and P are replaced by x and 
y 7 respectively. 

The function F(g) in Eq. I|25[) remains arbitrary, to 
first order in £ . It can be found from the analysis of 
the terms oc C~ 2 in the Lagrange equation for the op- 
timal path |l7j . which is cumbersome. The calculation 
can be simplified by noticing that F(g) determines the 
change of the quasienergy g averaged over vibration pe- 
riod 2ir/oj(g). From Eq. (|25|l . taking into account the 
explicit form of K^ ^ and expression l|21|l . we have 

(dg/dr) opt = -C l {xd x gW+yd y gW 



-F{g)u-\g) [(d x gWy- + (d y g^Y 
-K)" 1 Ug)M^(g) + F(g)N(°\g)) , (26) 



where overline denotes period averaging. 

Alternatively the evolution of quasienergy on the opti- 
mal path can be found using the Langcvin equation for 
dg/dr, which can be obtained by averaging Eq. over 
the period 2it/u){g). The resulting equation describes 
drift and diffusion of g. The optimal path for g can 
be obtained using the variational formulation similar to 
Eq. (O- Since the corresponding variational problem is 
onc-dimcnsional, it is easy to find the optimal path. As 
is often the case for one-dimensional systems driven by 



white noise, the optimal path is the time-reversed path 
in the absence of noise, 

(d#) opt = K)^,(p(»)(s). 

Comparing this expression with Eq. I)26|l we find 



F{g) = -2^{g)M^{g)/N^{g). 



(27) 



To lowest order in we can replace g with g on the 
optimal path. 

Equations (fTTft . (O, (|2SJ), and J2ZJ) describe the LS in 
a simple form of the Fourier transform of the functions 
F{g)d x ) ; x 1 F{g)d^y on the optimal path. The general ex- 
pression is further simplified near the peaks of %i{v). 
They occur where v is close to the frequency uj(g) or its 
overtones for certain values of g. It is convenient to write 
x and y as Fourier series, 



Ex(m; g) expfimip) 
71. 



and similarly for y. Then calculating the LS is reduced 
to taking the Fourier transform of the oscillating factors 
exp[im - op t(r)] weighted with functions of g pt(7") that 
smoothly depend on time. 

On the MPEP the leading term in the phase ip pt(T) 
is V'opt (t) ks J dr'tj(popt(T / )) . The function w(g) mono- 
tonically decreases with increasing g. Since gopt(j~) is 
monotonic as function of time, maj(g op t) can be in res- 
onance with v on the optimal path only at a certain 
time. This allows us to single out resonant contributions 
Xi(m; v) to Xi( u ) from the corresponding mth overtones 
of x{g,il)),y(g 1 ip). Near resonance integration over t in 
Eq. itT^ can be done by the steepest descent method, for 
slowly varying g(r). It gives 

Xi(m; v) = C m nC 1/2 \2uM^/\du/dg\ iV (0)2 l ^ 

x (y(-m; g um ), -x(-m; g^ m )j 

Here, C m is a phase factor, \C m \ = 1. The subscript g vm 
indicates that the expression in the brackets should be 
calculated for g = g vmi with g vm given by the condition 
mu){g mn ) = v. 

The change of the activation energy 8R\ is determined 
by the LS Xic{ v ) defined in Eq. (|17|) . From the explicit 
form of the matrix elements x(m; g),y(m; g) found in 
Ref.|9|it follows that arg[y(— m; g)x*(— m; g)} = —tt/2 for 
to > 0, in the considered range of g and (i. Then from 
the expression for Xi( m 'i l/ ) wc obtain that the spectral 
peak of xic near an mth overtone, xi c (m; v), is given by 
the expression 



1/2 



Xic(m; v) = irC 1/2 2vM^ /\du/dg\ A^ (0)s 
x(\x(-m;g vm )\ + \y(-m;g vm )\\ (to > 0). (28) 



The Fourier components \x(m; g)\, \y(m; g)\ rapidly de- 
crease with increasing |to| (exponentially, for large \m\ 
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0). Therefore the major peak of Xi c {v) is the peak of 
Xi c (l;^)- From Eq. 128fl . near the maximum it has the 
form 

Xi c (l;^)| oc (1 - az)0(z), z = w(g min ) - v, (29) 

where a ~ 1 is a numerical factor and 0(x) is the step 
function [u>(g m i n ) = 2(/i + l) 1 ' 2 ]. In obtaining this ex- 
pression we used that, for small g — ffmin, vibrations 
with given g are nearly sinusoidal, with x{\;g),y{\\g) oc 
(.9-.9min) 1/2 and M^ a \N^ oc g-g min - The sharp asym- 
metry of the peak of xi c (l; v) \ is due to the fact that the 
eigenfrequencies oj(g) have a cutoff at oj(g m i n )- 

A similar calculation shows that the peak from the 
2nd overtone is smoother, with Xi c (2;i/)| [2w(c/ m ; n ) — 
v] 1 / 2 for small 2o>(g m i n ) — v. The maximum of Xic(2; v) 
is shifted from 2oj(g m \- a ) to lower frequencies and has a 
smaller height than x\ c {l;v). Higher-order peaks have 
still smaller heights. 




FIG. 2: The multiple-peak LS Xic( J/ ) for an underdamped sys- 
tem. Solid lines show the overtones xi (m;f), Eq. J2HJ, with 
m = 1,2,3. Squares show the LS calculated by numerically 
finding the MPEP. The data refer to C = 5, fj, - = 1.5. 



The multi-peak structure of the LS for small damping 
is clearly seen in Fig. [2 The lines in this figure show 
asymptotic expressions (I28f) for the overtones X\ c {m\ v). 
The squares arc obtained by numerically finding the 
optimal path c^J t (t) and then calculating X\cW) fr° m 
Eqs. ifTTJl. |Q , and l(T7|l. As expected from Eq. (|2HJ), the 
major peak of the numerically calculated LS is at fre- 
quency w((? m i n ). Other peaks are located near overtones 
of w(g m in), on the low-frequency side. All peaks have 
characteristic strongly asymmetric shapes. 

Equation l|28|) for the major peak, m = 1, is in a good 
agreement with the numerical calculations. The agree- 
ment for the overtones is worse, because the peaks are 
much broader and the contributions to the LS from dif- 
ferent overtones overlap. The full LS is not given just 
by a sum over m of x\ c {m;v). It is necessary to take 
into account interference of the contributions from vibra- 
tions with different quasicnergies g vm but close mu)(g l/m ) 



and, of course, the effect of relaxation. We note that for 
v oc 2u>d — Lop < the peaks of the LS have much smaller 
amplitudes. 

The above results demonstrate that, for an under- 
damped oscillator, the rate of switching between period 
two states can be resonantly increased by applying an 
extra field with frequency u>d ~ u>f/2. The amplitude of 
the LS peaks is oc C _1 ^ 2 i it is parametrically larger than 
the LS at zero frequency, which determines the symmetry 
lifting and is oc C _1 - 

We note that, for extremely weak damping, a change 
of the switching rate with the field amplitude Ad may be 
associated with the field-induced mixing of the attraction 
basins, as in equilibrium systems (2(|. We do not discuss 
this mechanism here. 



VIII. CONCLUSIONS 

In this paper we considered an oscillator parametri- 
cally modulated by a comparatively strong field at nearly 
twice its eigenfrequency and additionally driven by a 
comparatively weak nearly resonant field. Because of the 
parametric modulation the oscillator displays period dou- 
bling. It has two vibrational states which differ only by 
phase in the absence of the additional field. Even a com- 
paratively weak additional field can strongly affect the 
oscillator by changing the rates of switching between the 
period two states. We have shown that the rate change 
depends exponentially on the ratio of the amplitude of 
this field Ad to the characteristic fluctuation intensity D. 
For small D this change becomes large even where the 
field only weakly perturbs the dynamics of the system. 

If the frequency of the additional field LOd coincides 
with the frequency of the period two states wf/2, the 
switching rates W12 and W21 between the states become 
different from each other. As a result, the stationary 
populations of the states become different, too. This is 
the effect of symmetry lifting. It depends exponentially 
strongly on the field amplitude. 

For small frequency difference, \u>d — Wf/2\ S W nm , 
the additional field leads to oscillations of the state 
populations which, in turn, lead to fluctuation-induced 
three- and multiple- wave mixing. For a larger detuning, 
\ujd — ^f/2| 3> W nm , the major effect of the additional 
field is the increase of the switching rates W nm averaged 
over the beat period 47r /\u>p — 2u>d\- 

Both for small and comparatively large \u>f — 2uid\ the 
change of the switching rates is characterized by the LS 
Xn- The latter gives the proportionality coefficient be- 
tween the field-induced change of the activation energy of 
switching from an nth state SR n and the field amplitude 
Ad- We have obtained an explicit general expression for 
the LS in terms of the path that the system is most likely 
to follow in switching, the MPEP. For ujd = ujf/2, the 
two components of the vector Xn give the coefficients in 
5R n at cos 4>d and sin^, where cf>d is the phase of the 
additional field relative to the phase of the strong field. 
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For \ujd — ojp/2\ 3> W nm the quantity of interest is \ic 
denned by Eq. (fT7|) . 

The major qualitative features of the LS are (i) scaling 
behavior near the bifurcation point where the period two 
states merge, and (ii) the occurrence of resonant peaks as 
a function of frequency for weak damping. We have found 
that x n scales with the distance rj to the bifurcation point 
a s Xn oc Thus the field-induced correction to the 

activation energy of switching decreases as the system 
approaches the bifurcation point. However, this decrease 
is significantly slower than that of the major term in the 
activation energy, which is oc rj 2 . 

Resonant peaks of the LS become pronounced where, 
in the rotating frame, oscillator vibrations about the pe- 
riod two states are undcrdamped. The vibration frequen- 
cies are much less than u>d, and the condition that these 
vibrations are undcrdamped is more restrictive than the 
requirement that the oscillator be undcrdamped in the 
laboratory frame. The LS displays several peaks. The 
major peak is shifted from cup/2 by the frequency of 
small-amplitude vibrations about the period two states. 
Other peaks are shifted approximately by the overtones 
of this frequency and have smaller amplitudes. All LS 
peaks have characteristic strongly asymmetric shapes. 

The effects discussed in this paper are not limited to 
a parametric oscillator, they can be observed in other 
systems with period two states. The results on mod- 



ulation of switching rates by an additional field can 
be extended also to other systems with coexisting vi- 
brational states, including resonantly driven nonlinear 
nano- and micromechanical oscillators and Josephson 
junctions. Fluctuational interstate transitions in these 
systems were recently studied experimentally |2^ . l28ll29j . 
Although the theory of the LS for modulated oscillators 
was not developed until this paper, a resonant increase 
of the switching rate by an additional field was expected 
from qualitative arguments based on the analogy with 
static systems |30j . The preliminary experimental data 
indicate that the effect occurs in Josephson junctions 
when the frequencies of the both strong and weak fields 
are close to the plasma frequency |U . 

In conclusion, we have studied the effect of an addi- 
tional field on a fluctuating paramctrically modulated os- 
cillator. We predict strong change of the populations of 
the period two states by a comparatively weak field. We 
also predict that the logarithm of the rate of interstate 
switching is linear in the field amplitude, and the pro- 
portionality coefficient may display resonant peaks as a 
function of the field frequency 
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